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1.1 Introduction: The high power response of the transmon-cavity 
system 

The Jaynes-Cummings (JC) Hamiltonian provides a quantum model for a two-level 
system (qubit) interacting with a quantized electromagnetic mode. It is widely appli- 



cable to experiments with natural atoms (Haroche and Raimond, 2006 


Gleyzes et a/., 


2007| Boca et a/., 2004 


Brennecke et a/., 2007 Maunz et al, 2005) as well as for solid- 


state 'artificial atoms' (Reithmaier et a/., 2004 


Yoshie et a/., 2004 Wallraff et a/., 



2004). The discussion in this chapter applies to any realization of the model that can 
reach the appropriate parameter regimes and be driven sufficiently strongly, but for 
concreteness we adopt the language and focus on typical parameters from the field of 
circuit quantum electrodynamics (circuit QED), where the relevant parameter range 
is easily achieved in experiments. 

We write the Jaynes-Cummings model with drive and dissipation (ft = 1) 



H = u c a^a - 



g(a + a 



ffl ) + H~f + H K 



(1.1) 



with cavity frequency cj c /27t, qubit frequency Co> q /27r and where is the time- 
dependent drive of the cavity, g is the cavity-qubit coupling, and H liK represent the 
coupling to the qubit and cavity baths, respectively. 

The JC Hamiltonian can be diagonalized analytically, but in the presence of a 
drive and dissipation the open-system model becomes non-trivial, with the ef- 
fective behavior depending strongly on the specific parameter regime. Several central 
parameters are involved in the classification of theses different regimes. The case where 
the cavity relaxation rate k is much larger (smaller) than the two-level dissipation and 
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dephasing rates 7, 7^ is known as the bad (good) cavity limit. The strong dispersive 



regime (Gambetta et a/., 2006 Schuster et a/., 2007) of the JC model describes the 



situation that the presence of the qubit causes the cavity frequency to be shifted by an 
amount x much greater than the cavity linewidth. The shift depends on the number 
of excitations in the systems x = x(N) and for low photon numbers is similar to a 
Duffing oscillator nonlinearity . 

In cavity QED, various forms of "single atom bistability" in this model are known: 



single atom absorptive bistability (Savage and Carmichael, 19881 Kerckhoff et 



2011), exists in the weak coupling regime in the good cavity limit < 7 < ft); also 



closely related is the single atom phase bistability of spontaneous dressed state po- 
larization (Alsing and Carmichael, 1991 Kilin and Krinitskaya, 1991] ) which concerns 
the case where the atom and the cavity are in resonance 5 = ou c — ou q = 0. 

In this chapter we analyze the high excitation nonlinear response of the Jaynes- 
Cummings model in quantum optics when the qubit and cavity are strongly coupled. 
We focus on the parameter ranges appropriate for transmon qubits in the circuit quan- 
tum electrodynamics architecture, where the system behaves essentially as a nonlinear 
quantum oscillator and we analyze the quantum and semi-classical dynamics. One of 
the central motivations is that under strong excitation tones, the nonlinear response 
can lead to qubit quantum state discrimination and below we present initial results for 
the cases when the qubit and cavity are on resonance or far off-resonance (dispersive). 



1.2 Implications of the nonlinearity at the high excitation regime 

A characteristic feature of the JC model is that for very high excitation number N ^> 
1, the excitation number-dependent frequency shift obeys x(iV) — »• 0: the transition 
frequency returns to the bare cavity frequency. In the presence of dissipation this 
happens effectively when x(iV) < ft, and for all larger N the response of the system 
is linear with respect to the drive. We describe this behavior as setting in at an 
excitation number Abare? with the definition x(^bare) = ft. In the strong dispersive 
regime we have A^are ^> iV crit , where A^t as usual denotes the excitation level where 
the dispersive approximation breaks down (defined below). The latter inequality has an 
important consequence for the theory: a perturbative expansion in the small parameter 



A^/A^rit, typically useful ( Boissonneault et a/., 2008 Boissonneault et a/., 2009) in the 



dispersive regime, is not applicable for the interesting regime N > Abare where the 
system regains the linear response. 

We discuss two different parameter regimes, the semi-classical and quantum de- 
generate, shown schematically in fig. |1.1| We note that a very related semi-classical 
nonlinearity can be introduced via adding a Josephson junction directly to the res- 
onator. This approach leads to the Josephson Bifurcation Amplifier (JBA), the Cavity 
Bifurcation Amplifier (CBA) which realize an efficient qubit readout based on the 
dynamical bistability of the resonator state, and were the first to make use of bright 
states for qubit measurement (|Siddiqi et a/., 2006} |Boulant et a/., 2007 Metcalfe et a/., 



2007). In this review we confine ourselves to the minimal model of only one qubit in- 
teracting with one linear cavity which has a different structure of anharmonicity with 
unique effects. 




Fig. 1.1 Qubit state measurement in circuit QED can operate in different parameter regimes 
and relies on different dynamical phenomena of the strongly coupled transmon-resonator 
system. The dispersive readout is the least disruptive to the qubit state and it is realized 
where the cavity and qubit are strongly detuned. The high power readout operates in a 
regime where the system response can be described using a semi-classical model and yields 
an relatively high fidelity with simple measurement protocol. When the cavity and qubit are 
on resonance (the quantum degenerate regime) it is theoretically predicted that the photon 
blockade can also be used to realize a high fidelity readout. 



In Section 1.2.1} we consider the model in the bad cavity limit and on timescales 
short compared to the atomic coherence time where the dynamics are those of a nonlin- 
ear oscillator, which we study by both semiclassical methods and quantum trajectory 
simulations. Our main result is that there exists a threshold drive £02 at which the 
photon occupation increases by several orders of magnitude over a small range of 
the drive amplitude. We perform both quantum trajectory simulations and a non- 
perturbative semiclassical analysis, including the drive and the cavity damping. Our 
results are in qualitative agreement with recent experiments (Reed et a/., 20101) for 



a circuit quantum electrodynamics (QED) device (DiCarlo et a/., 2010) containing 4 



transmon (Koch et a/., 2007 Schreier et a/., 2008|) qubits, demonstrating that the JC 



model captures the essential physics despite making an enormous simplification of the 
full system Hamiltonian. 



In Section 1.2.2 the situation where the qubit and the cavity are on resonance is 



shown to lead to a coexistence of photon blockaded states and highly excited quasi- 
coherent states (QCS) with the same driving tone. This is also a result of the nonlin- 
earity that arises in the Jaynes-Cummings model, but in contrast this regime cannot 
be fully described with a semiclassical theory and requires exact quantum simulations 
to reveal. 

1.2.1 Transient response in the dispersive regime 

The behavior of the JC nonlinearity goes beyond the Kerr nonlinearity that is often 
considered. Dispersive bistability ( jMarburger and Felber, 1978| ) from a Kerr nonlinear- 
ity is well-known in atomic cavity QED (Gibbs et a/., 1976| ). It has been implemented 



in the solid state via the nonlinearity of a Josephson junction (Siddiqi et a/., 2004), 



and in the circuit QED architecture has produced high-fidelity readout of qubits ( Sid- 



diqi et a/., 2006 Boulant et a/., 2007 Mallet et a/., 2009). Similar schemes have been 
implemented with nonlinear micromechanical resonators (Almog et a/., 2006[). How- 



ever, unlike the Kerr anharmonicity, the JC anharmonicity does not remain constant 
but rather diminishes toward zero as the cavity occupation is increased. As a result, 
for sufficiently strong drive the response of the JC model must return to the linear 
response of the bare cavity. Instead of coherent driving, an alternative way to saturate 
the qubit and cause the JC system response to return to the bare cavity response 
is to couple the system to a bath at elevated temperature, as has been investigated 



theoretically (Dykman and Ivanov, 1976 Rau et a/., 2004) and experimentally (Fink 



et al, 2010P - 

Recent experiments ( Reed et a/., 2010[ ), which operate in both the strong dispersive 
regime and the bad cavity limit, show a nontrivial response under conditions of strong 
drive, arousing interest due to its usefulness for high-fidelity qubit readout. 

For driving at a single frequency, we write the driven JC Hamiltonian 

H = uj c o)a + ^(Jz + + aV_) + -4= (a + a)) cos(cjd^), 

2 v2 

given with the drive frequency uo&. Operating in the strong-dispersive bad-cavity regime 
defines a hierarchy of scales 

7, 70 < « < g 2 /$ < # < £ < u; c , (1.2) 

where 5 = uj q — uj c is the qubit- cavity detuning. We can make the standard transfor- 
mation (Carbonaro et a/., 1979) H = T~ 1 HT to decouple the qubit and cavity 



H = uj c a^a + (uj c — A)- 



(a + a T ) cos(co>d£) 5 



(1.3) 



2 ' V2 

dropping terms from the transformed drive that are suppressed as 0{N~ 1 / 2 ) and 
0{g/S). The resulting Hamiltonian would be trivial were it not for the fact that the 
transformation T, defined by 

T = exp[-6>(4A0- 1/2 (a<7 + + a* a J)], (1.4) 

sin 9 = -2 5 iV 1 /2/A, cos 9 = 6/ A, (1.5) 

A = (S 2 +4g 2 N) 1/2 , (1.6) 
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Fig. 1.2 Transmitted heterodyne amplitude | (a) | as a function of drive detuning (normalized 
by the dispersive shift % = g 2 /5) and drive amplitude (normalized by the amplitude to put 
n — 1 photon in the cavity in linear response, £i = K,/y/2). Dark colors indicate larger 
amplitudes, (a) Experimental data (Ree d and DiCarlo, 2010| , for a device with cavity at 
9.07 GHz and 4 transmon qubits at 7.0,7.5,8.0,12.3 GHz. All qubits are initialized in their 
ground state, and the signal is integrated for the first 400 ns ~ 4/k after switching on the 
drive, (b) Numerical results for the JC model of eqn |1.7| with qubit fixed to the ground state 
and effective parameters 6/2tt = -1.0 GHz, g/2n = 0.2 GHz, k/2tt = 0.001 GHz. These are 
only intended as representative numbers for circuit QED and were not optimized against the 
data of panel (a). Hilbert space is truncated at 10,000 excitations (truncation artifacts are 
visible for the strongest drive), and results are shown for time t = 2.5/ k. 

depends on the total number of excitations, N = a^a + a z /2 + 1/2. For photon decay 
at rate n we can write the decoupled quantum master equation after dropping small 
terms, 

p = -i[H, p] + K ([ap, at] + [a, pa)))/2, (1.7) 

which we integrate numerically in a truncated Hilbert space using the method of quan- 
tum trajectories, after making the rotating wave approximation (RWA) with respect 
to the drive. The experiments we wish to describe are performed on a timescale short 
compared to the qubit decoherence times 7 -1 , and we therefore treat a z as a 
constant of motion. The remaining degree of freedom constitutes a Jaynes-Cummings 
oscillator. Note that the qubit relaxation and dephasing terms that we have dropped 
involve the <j± and <j z operators and would transform in a nontrivial way under the 



decoupling transformation T ( Boissonneault et a/., 2008), ( Boissonneault et a/., 2009) 
The results of the numerical integration for <j z = — 1 are compared with recent ex- 



perimental data (Reed and DiCarlo, 2010) in Fig. 1.2, where we show the average 



heterodyne amplitude \{a)\ as a function of drive frequency and amplitude. Despite 
the presence of 4 qubits in the device, the fact that extensions beyond a two- level 
model would seem necessary since higher levels of the transmons are certainly occu- 
pied for such strong driving] and despite the fact that the Rabi Hamiltonian might 



1 Simulations show approximately 10 transmon levels are required to simulate the experiment 
quantitatively. 



seem more appropriate for such large photon occupation, y Ng ~ cj c , nevertheless the 
JC model qualitatively reproduces the features of the experiment^ In particular, for 
weak driving we see a response as expected at the dispersively shifted cavity frequency 
uj c — x, with x — 9 2 A 5 which shifts towards lower frequencies as the drive increases. For 
stronger driving a dip appears in the response, which we interpret as a consequence of 
plotting the absolute value of the ensemble- averaged amplitude (a) in the classically 
bistable region, as we discuss below. For increasing drive the dip shifts to lower fre- 
quencies; finally for the strongest driving, the response becomes centered at the bare 
cavity frequency uj c /2tt and is single-peaked and extremely strong. We note that both 
the experiment and numerical integration are terminated at a transient time of only a 
few cavity lifetimes, and we have checked that the numerical response is substantially 
different for the steady state (see below). 

When there is a large number of photons in the system such that the anharmonicity 
is greatly diminished, it is possible to use a semiclassi cal model, similar to (|Alsing 
and Carmichael, 1991), (Kilin and Krinitskaya, 1991), (Peano and Thorwart, 2010), 



to characterize the transmission. In fact this is also a good approximation for the 
response at low powers in the dispersive regime when the ratio of the anharmonicity 
of the dispersive Hamiltonian to the decay rate (width of the levels) is such that 
the N — 1 -<->_7V photon peak overlaps well with the N «->• N + 1 photon peak. By 
to second order in g 2 N/S 2 this condition can be seen to be N ^> iV sc , 
(for the parameters of Fig. 



expanding 
where 7V SC 



1.3 



1.2 3, N sc = 1.6). In the opposite limit 



we will see photon blockade and associated effects, as in (Bishop et a/., 2009|). Recently 



it was shown that it is possible to have a coexistence of both the semiclassical and 



quantum solutions for a certain range of parameters of the system and drive ( Ginossar 
\et a/., 2010) as we discuss below in Section|1.2.2l The semiclassical model will remain 



valid for N > iV cr it, where a perturbative expansion of the Hamiltonian 1.3 in term s of 
N/N cr [ t fails to converge, where iV cr it = 5 2 /Ag 2 . We rewrite the Hamiltonian eqn 
using canonical variables X — y/l/2(a) + a) and P = \^/l/2(a) — a), 



1.3 



UJr 



H=^(X 2 



a z ) + ^X cos(u d t) 



^-^2g 2 (X 2 +P 2 +a z )+5 2 . 



(1.8) 



The semiclassical approximation consists of treating X and P as numbers. Instead of 



directly solving Eq. 1.7 we write the equations of motion for X, P from the diagonalized 
Hamiltonian H of the closed system and the effect of cavity relaxation is incorporated 
through a phenomenological damping term proportional to k. We solve for the steady 
state, treating X 2 + P 2 as a constant (thus we ignore harmonic generation), giving a 
nonlinear equation for the amplitude A = \/X 2 + P 2 



2 c2 



A 2 = 



(1.9) 



2 We emphasize that the effective parameters in the simulation are of the same magnitude as in 
the experiments but we do not expect any quantitative correspondence. 



with amplitude-dependent frequency shift x(-A), given b}J^] 

m -'- wmm* - <L10) 

This reproduces for small driving the usual dispersive shift x(0) ~ ±# 2 /£ and for 
strong driving shows the saturation effect lim^^oo x(A) = 0. 

The solution of eqn L9 is plotted in Fig. L3 for the same parameters as in Fig. |1.2] o. 



For weak driving the system response approaches the linear response of the dispersively 
shifted cavity. Above the lower critical amplitude £ci the frequency response bifur- 
cates, and the JC oscillator enters a region of bistability. We denote by C\ the point at 
which the bifur catio n first appears. Dropping terms which are small according to the 



hierarchy of eqn 1.2 this point occurs at £ci = (£ft) 3 / 2 3 -3 / 4 g -2 , Vtci — x(0) — a/3^/2, 
writing the drive detuning as Q = uod — uo c . The dip in the heterodyne measurement 
of Fig. 1.2 appears within the bifurcation region (Fig. |1.3fr ), indicating that this dip 



is the result of ensemble- averaging of the coherent heterodyne amplitude in the region 
of classical bistability. In Fig. |1 .3| o we see that the semiclassical and quantum simu- 
lation yield the same response outside the region of bistability. Within the region of 
bistability, quantum noise causes switching ( |Dykman and Smelyanskiy, 1 988) between 
the two semiclassical solutions, one dim and one bright, with almost opposite phases. 
An analytical derivation of the dip in the steady-state amplitude for a Kerr nonlin- 



earity was given in (Drummond and Walls, 1980|). In our case, both the experiment 



and numerical integration are terminated at a transient time of only a few cavity life- 
times. Therefore the exact form of the averaged response is influenced by the initial 
conditions. We have checked that the position of the dip in the numerical response 
is shifted towards lower frequencies in the steady state, consistent with the switching 
rate being slow compared to the cavity decay rate. 

As the drive increases, and unlike the Kerr oscillator, the frequency extent of 
the bistable region shrinks and eventually vanishes at the upper critical amplitude 
ic2 — g/V%- In the effective theory the upper critical point C2 is located very close to 
the bare cavity frequency. This indicates that for driving at the bare cavity frequency, 
there is no bistability, but rather a finite region (a 'step') in the vicinity of the critical 



point (Fig. 1.3 3), where the response becomes strongly sensitive to the drive amplitude. 
The size of the step can be shown to be a factor of r = Abright/Aiim = 2g 2 / n5 in 
amplitude, and represents a very high gain (dA/d£ = y/2g/ ft 3//2 (5 1//2 ) in the strong- 
dispersive regime. Above the step we see that the response approaches the linear 
response of the bare cavity as N ~ A^are- 

1.2.2 The quantum degenerate regime 

We now move to exploring what happens when the detuning between the qubit and 
the cavity is reduced such that the anharmonicity of quantum ladder of states becomes 
much larger than the corresponding linewidth n (see Fig. |1.1[ ). In order to describe 
the response of the system to external drive in this regime it is important to take into 
account the quantum dynamics on the lower anharmonic part of the ladder. When the 



3 Note that TV w A 2 in the semiclassical approximation and in that case x(iV) = x(A 2 ) = x(A). 




Fig. 1.3 Solution to the semiclassical equation |1.9| using the same parameters as Fig. |1.2| 3. 
(a) Amplitude response as a function of drive frequency and amplitude. The region of bifurca- 
tion is indicated by the shaded area, and has corners at the critical points Ci, Ci. The dashed 
lines indicate the boundaries of the bistable region for a Kerr oscillator (Duffing oscillator), 
constructed by making the power-series expansion of the Hamiltonian to second order in 
N/N cr it. The Kerr bistability region matches the JC region in the vicinity of C\ but does not 
exhibit a second critical point, (b) Cut through (a) for a drive of 6.3£i, showing the frequency 
dependence of the classical solutions (solid line). For comparison, the response from the full 
quantum simulation of Fig. |1 .2| o is also plotted (dashed line) for the same parameters, (c) Cut 
through (a) for driving at the bare cavity frequency, showing the large gain available close 
to C2 (the 'step'). Faint lines indicate linear response, (d) Same as (c), showing intracavity 
amplitude on a linear scale. 



system is initialized in the ground state, there is a range of drive strengths for which the 
system will remain blockaded from excitations out of the ground state. However, since 
the anharmonicity of the JC ladder decreases with excitation number, the transition 
frequency for excitations between adjacent levels ultimately approaches the bare cavity 
frequency. Qualitatively, when the excitation level n is such that the anharmonicity 
becomes smaller than the linewidth ft, we expect the state dynamics to be semiclassical, 



similar to a driven-damped harmonic oscillator (Alsing and Carmichael, 1991 



Kilin 



and Krinitskaya, 1991). In this regime we therefore expect the drive to excite states 
which have an occupation function similar to coherent states but are highly mixed, 
which we call quasi-coherent states (QCS). More specifically, in order to support a 
coherent wave packet centered around level n, with a standard deviation of y/n, the 
difference of transition frequencies across the wavepacket has to be of the order of the 
linewidth ft. This approximate criterion for a minimal n can be written as co> n +2cr — 
Un-2a ~ k where uj n ±2 a are the ladder transition frequencies, positioned 2a = 
above and below the mean level n. 

To study this system, we use the stochastic Schrodinger equation which describes 
the quantum evolution of the qubit-cavity system in the presence of qubit and cavity 
dissipation and a direct photodetection model ( |Breuer and Petruccione, 2QQ2[ ) 

dip(t) = -i (H(t) - ~ \ a ^ a -) ^)dt + (l.H) 

a - m m)dN 1 + (^K-m)dN 2 (i.i2) 



with the H(t) given by Eqn. 1.2 Using numerical integration of this equation we 
generate an ensembles of trajectories of the quantum evolution of the wave function 
conditioned on the measurement signal. Quantitatively, we find that the lifetime of 
the QCS is long but finite, and increases with the amplitude of the drive (see inset in 



Fig. 1.4). 



As we explain below, we find that low-lying QCS (n = 20) are the most effective for 
optimizing the overall readout fidelity. Note that the JC ladder consists of two mani- 
folds (originating from the degeneracy of the bare states \g, n + 1) and |e, n)) denoted 
by (±), and we will always refer to states occupying one manifold since the drive is off- 
resonant with respect to the other manifold. Transitions between manifolds contribute 
to the decay of the QCS to the dim state. Such transitions can be induced by the drive 
but their rate is smaller by a factor of (9(n -1 / 2 ) compared to the rate of transitions 
inside the same manifold. An additional source of inter-manifold transitions are decay 
(Ti) and pure qubit dephasing (T^), whose effects in the presence of drive were studied 
in the context of the dispersive regime (|Boissonneault et a/., 2009). Indeed, as we see 



in Fig. |1.4[ changing T x has a noticeable effect on the QCS lifetimes. For very large n, 
these processes become ineffective for inducing decay of QCS, since then the difference 
between the manifold excitation frequencies becomes smaller than ft, and therefore the 
drive effectively drives both manifolds. For superconducting transmon qubits 7\ is the 



dominant decay process, and we show its effect on the overall fidelity in Fig. [L8 

The QCS exist with drive amplitudes where the ground state is photon blockaded, 
giving rise to a dynamical bistability between quantum and semi-classical parts of 
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Fig. 1.4 Schematic figure showing an interpretation of the numerical simulations. Long-lived 
quasi-coherent states (QCS) decay due to photon emission events which occur in the expo- 
nential tail of the wave packet in the lower anharmonic parts of the ladder of states of the 
Jaynes-Cummings model. These emissions are uncompensated by the detuned drive and as 
a result the wave packet eventually falls into the blockaded regime (n) —> 0. The inset shows 
the results of quantum simulations: the lifetimes (r&) of QCS is increasing with drive strength 
and mean photon number for drive strengths where the dim state is still photon blockaded 
demonstrating the coexistence regime. We also see that the qubit decay (Ti) has a distinct 
influence on the lifetime of the QCS (here the coupling was taken to be g/27i — 100MHz and 
the cavity lifetime here is ft -1 = 60ns). 

the JC ladder. Indeed, we see that there is a basin of attraction for states initialized 
as coherent wave packets to persist as QCS, and we characterize it according to the 
probability of the state to decay on the timescale 75 ^> n~ x . It is worth noting that 
this metastability cannot be described using a semiclassical picture in cases where a 
significant part of the Hilbert space where the tail of the QCS resides is quantized 
(i.e. "low lying" QCSs). Similarly processes with the opposite transition where the 
system succeeds in leaving of the blockade into the QCS basin cannot be described 
without solving the full quantum dynamics of the wave packet. The characterization 
and exact mechanisms of dynamical transitions between these states in this regime 
are a matter of ongoing research. Recent experiment have detected exponentially long 
lifetimes for quasi- coherent states excited by a chirp in a high excitation state of a 



Josephson phase qubit (Shalibo et a/., 2011) and driven by a holding tone. In Fig. 1.5 



we plot the contours of equal probability of the QCS to decay to a manifold of states 
close to the ground states, after a time ft -1 , given that it was initialized with a certain 
amplitude (a) and phase (0). We see a large region supporting QCS, and the phase 




Initial state amplitude a 

Fig. 1.5 Probability for the QCS to decay after being initialized as coherent state wave 
packet \ae l9 ) and driven for a time k~ x . The equal probability contours trace out two basins 
of attraction: states initialized inside the low decay probability contour (Pd < 1%) end up 
long lived (r& ^> ft -1 ), whereas states initialized left of Pd — 90% quickly decay to the ground 
state and remain photon blockaded. The parameters are g/2n — 100MHz, k/2ix — 4.05MHz, 
£/2tt = 9.9MHz,(cj d - cj c )/2tt = 12.3MHz,(cj c - u q )/2n = 9.7MHz, Ti = 1591ns (uo d is the 
drive frequency). 

sensitivity can be understood qualitatively from the time dependent simulations: a 
mismatch between the phases of the drive and initial coherent state causes ringing of 
the wave packet outside of the basin into the too anharmonic part of the ladder from 
which it cannot recover. In addition to the existence of this basin, the anharmonicity 
acts together with the cavity decay to induce mixing of the QCS: even for bright states 
(n > 40) we extracted a relatively low purity of Tr(p 2 ) < 0.5. 

1.3 Applications for high fidelity qubit state measurement 

The nonlinear response either in the dispersive or degenerate regimes opens the possi- 
bility of using the system's own high gain for (self) qubit readout. The figures of merit 
of a useful scheme involve low error rate, high contrast, speed and ideally the property 
of quantum non-demolition. To these we should probably add robustness because for a 
protocol to be useful it should sustain experimental imperfections. Fundamentally, the 
source of the high gain in the two regimes is in the nonlinearity of the two-level system 
(usually the two lowest levels of a quantized anharmonic oscillator). The meaning of 
nonlinearity and high gain for the semiclassical vs. quantum degenerate regimes are 
quite different and will be discussed below. What is common in these two schemes is 
that both require the application of relatively strong tones, they both make use of a 
large part of the phase space of the system, and both are projective. In addition both 
schemes have many degrees of freedom in frequency and amplitude modulation of their 



control pulses and are therefore amenable to optimization. For the scheme operating in 
the quantum degenerate regime we have implemented a partial optimization in order 
to demonstrate that high fidelity readout is feasible. 



1.3.1 Symmetry breaking for the transmon device 



From the semiclassical eqns |1.9[ |1 . 10 it follows that for A ^> 1 the response of the 



system will have an approximate symmetry of reflection with respect to the bare cav- 
ity frequency A(Q,a z = +1) « A(— fl,cr z = —1). Therefore the response at the bare 
cavity frequency will be nearly independent of the state of the qubit, with respect to 
both the low and high power regimes. In order to translate the high gain available at 
the step into a qubit readout, it is necessary to break the symmetry of the response of 
the system between the qubit ground and excited states, such that the upper critical 
power £c2 will be qubit state dependent. In the JC model the symmetry follows from 
the weak dependence of the decoupled Hamiltonian H on the qubit state for high 
photon occupation. However, the experimentally-observed state dependence may be 
explained by a symmetry breaking caused by the higher levels of the weakly anhar- 
monic transmon, or by the presence of more than one qubit, see Fig. |1.6| Such an 
asymmetric qubit dependent response is well known for the transmon in the low power 



dispersive regime (Koch et a/., 2007). For the high power regime, such a dependence 



2010) by taking into account the full nonlinearity 



was shown in ( jBoissonneault et 
of the transmon. 

When designing a readout scheme that employs such a diminishing anharmonic- 
ity, the contrast of the readout is a product of both the symmetry breaking and the 
characteristic nonlinear response of the system near the critical point Ci. Experi- 
ments ( |Reed et a/., 20 10) initially were able to use this operating point to provide 
a scheme for qubit readout, which is attractive both because of the high fidelities 
achieved (approaching 90%, significantly better than is typical for linear dispersive 
readout in circuit QED ( jWallraff et a/., 200 5; Steffe n~et a/., 20ld| ) and because it does 
not require any auxiliary circuit elements in addition to the cavity and the qubit. This 



nonlinear response was also used for characterizing three-qubit GHZ states (DiCarlo 



et a/., 2010) and very recently for the characterization of transmon qubits designed in 
a three dimensional architecture ( Paik et a/., 201 1) . 



1.3.2 Coherent control in the quantum degenerate regime 

Qubit readout in solid state systems is an open problem, which is currently the subject 
of intensive experimental and theoretical research. High-fidelity single-shot readout is 
an important component for the successful implementation of quantum information 



protocols, such as measurement based error correction codes (Nielsen and Chuang, 



2000 ) as well as for closing the measurement loophole in Bell tests ( Garg and Mermin, 



1987 



Kofman and Korotkov, 2008 Ansmann et al, 2009). For measurements where 



the observed pointer state depends linearly on the qubit state, for example dispersive 



readout in circuit QED (cQED) (Blais et a/., 2004), there exists a unified theoretical 



understanding (Clerk et a/., 2010). Experimentally, these schemes require a following 
amplifier of high gain and low noise, spurring the development of quantum limited 
amplifiers (Spietz et a/., 2009 Bergeal et a/., 2010; Castellanos-Beltran et a/., 2008). 




Fig. 
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Symmetry breaking. State-dependent transition frequency 
En, q ) versus photon number n, where E n , q denotes energy of the 



system eigenstate with n photons and qubit state q: (a) for the JC model, parameters 
and |1.3| (b) for the model extended to 2 qubits, 5i/2tt = —1.0 GHz, 
6/2/27T = 0.25 GHz. Here, \2 denotes the 0-photon dispersive 



1.2 



as in Figs 

S 2 /2tt = -2.0 GHz, gi /27T 
shift of the second qubit; (c) for the model extended to one transmon qubit (Koch et al, 
2007|, tuned below the cavity, cj c /2tt = 7 GHz E c /2n = 0.2 GHz, Ej/2tt = 30 GHz, 

-0.5 GHz, 6 12 /2tv = -0.7 GHz, 



7 GHz E c /2tt 

g/2n = 0.29 GHz. (For the given parameters, 5oi/2ty = 



defining Sij = Ej — Ei 



with Ei the energy of the ith transmon level.) In all panels, 



the transition frequency asymptotically returns to the bare cavity frequency. In (a) the 
frequencies within the a z — ±1 manifolds are (nearly) symmetric with respect to the bare 
cavity frequency. For (b), if the state of one ('spectator') qubit is held constant, then 
the frequencies are asymmetric with respect to flipping the other ('active') qubit. In (c), 
the symmetry is also broken due the existence of higher levels in the weakly anharmonic 
transmon. 

However, the highest demonstrated fidelities to date rely on nonlinear measurement 
schemes with qubit dependent latching into a clearly distinguishable state, e.g. Joseph- 



son Bifurcation Amplifier ( JB A) as well as optimized readout of phase qubits ( Siddiqi 



et ah, 2006| |Mallet et a/., 2009| |Ansmann et a/., 2009[ ). For this class, during the mea- 
surement, the system evolves under the influence of time varying external fields and 
nonlinear dynamics, ultimately projecting the qubit state. The space for design and 



control parameters is very large, and the dependence of the readout fidelity on them is 
highly nontrivial. Therefore the optimization is difficult and does not posses a generic 
structure. 

We propose a coherent control based approach to the readout of a qubit that 
is strongly coupled to a cavity, based on an existing cQED architecture, but not 
necessarily limited to it. This approach is in the spirit of the latching readout schemes, 
but it differs in that the source of the nonlinearity is the Jaynes-Cummings (JC) 
interaction. When the qubit is brought into resonance with the cavity mode, the strong 
anharmonicity of the JC ladder of dressed states can prevent the excitation of the 
system even in the presence of a strong drive, a quantum phenomenon known as photon 
blockade (Imamoglu et a/., 1997 Birnbaum et a/., 2005 Hoffman et a/., 2011| ). However, 
due to fact that the JC anharmonicity is diminishing with the excitation number, we 
find a form of bistability, where highly excited quasi-coherent states (QCS) co-exist 
with the blockaded dim states (Fig. 1.5). In order to make use of this bistability to 
read out the qubit, it is necessary to solve the coherent control problem of selective 
population transfer, which is how to steer the system towards either the dim state 
or the QCS, depending on the initial state of the qubit (Fig. |1.7[ ). This selective 
dynamical mapping of the qubit state to the dim/bright states constitutes the readout 
scheme. It is potentially of high contrast, and hence robust against external amplifier 
imperfections. An advantage of this readout is that it uses no additional components 
beyond the qubit and the cavity, both already present as part of the cQED architecture. 
Based on a full quantum simulation which includes dissipation of the qubit and the 
cavity (we ignore pure dephasin^), we predict that implementing this scheme should 
yield very high fidelities between 90% and 98% for a typical range of realistic cQED 
parameters (Fig. |1.8| ). 

In the regime of coexistence, the dim quantum state and QCS present us with the 
possibility of implementing a high contrast readout scheme. This requires the solution 
of the coherent control problem of steering the logical 1 1) state to some point on 
within the basin of attraction (Fig. 1.5), while keeping the |^) far from the basin, in 
the manifold of dim states. In the presence of dissipation, the latter would quickly 
decay to the ground state, and remain there even in the presence of driving, due 
to the photon-blockade, whereas the QCS would persist for a long time T5 and emit 
approximately ft (71)575 photons. The standard coherent control problem of population 
transfe r (|Bergmann et a/., 1998 ), which was also discussed recently for superconducting 
qubits ( |Jirari et a/., 2009 ), is to maximize the probability P%^f of steering the state 
\i) to the state |/). However, here the goal is to bring the probability for selective 
steering Pi^f + Pi'^f < 2 close to its theoretical maximum, which is an essentially 
different coherent control problem. For systems with very large anharmonicities, for 
example atomic systems it is possible to effectively implement a population transfer 
via adiabatic control schemes such as STIRAP ( jBergmann et a/., 1998 ). The JC ladder 
anharmonicity is relatively small compared to atomic systems, and so these schemes 
are inapplicable here. 



4 Typically transmon qubits operate in a regime of Ej/Ec ^> 1 where pure dephasing is exponen- 
tially suppressed. 




Fig. 1.7 (Color) Readout control pulse (a) Time trace of the drive amplitude: a fast initial 
chirp (10 ns) can selectively steer the initial state, while the qubit is detuned from the cavity 
((cj q —cj c )/27t ~ 2g). It is followed by a slow displacement to increase contrast and lifetime of 
the latching state, while the qubit is resonant with the cavity (k/2tt — 2.5 MHz). The drive 
amplitude ramp is limited so that the photon blockade is not broken, but the contrast is 
enhanced by additional driving at the highest drive amplitude, (b) A diagram of transition 
frequencies shows how the drive frequency chirps through the JC ladder frequencies of the (+) 
manifold, and how the manifold changes due to the time dependent qubit frequency, (c) Wave 
packet snapshots at selected times (indicated by bullet points on panel (b)) of the chirping 
drive frequency of panel (b) conditioned on the initial state of the qubit. (d) The temporal 
evolution of the reduced density matrix |p mn | (the x,y axes denote the quantum numbers 
m, n of the cavity levels) of the cavity with the control pulse (a) when the qubit initial 
state is superposition -^=(|0) + |1)). The resonator enters a mesoscopic state of superposition 
around t = t c due to the entanglement with the qubit and the quantum state sensitivity 
of the protocol. At later times the off-diagonal parts of this superposition dephase quickly 
due to the interaction with the environment and the state of the system is being completely 
projected around t — 3t c - 



The control pulse sequence we apply is depicted in Fig. |1.7[ and consists of three 
parts: (1) a strong chirped pulse (t < t c ) drives the cavity and passes through the 
resonance of the cavity, with the qubit being detuned. Due to the interaction with the 
qubit, the cavity behaves as nonlinear oscillator with its set of transition frequencies 
depending on the state of the qubit (see the two distinct sets of lines in Fig. |1.7[ b)). 
The cavity responds with a ringing behavior which is different for the two cases (see 
Fig. |1.7[ c)). The ringing due to the pulse effectively maps the \\) and |t) to the dim and 
bright state basins, respectively (see Fig. 3(c)). Since Kt c <C 1, an initial superposition 
a |t) + f3\l) maps into a coherent superposition of the dim and bright states. Next, 
(2) a much weaker long pulse transfers the initially created bright state (for initial 
|t)) to even brighter and longer lived states (t c < t < t^) , and (3) steady driving for 
additional contrast (th < t < tf). For an initial superposition the interference terms 
between the dim and bright states decohere on the timescale of ft -1 , such that the 
interaction with the reservoir for t > t c effects a projection of the pointer state. In 
designing such a pulse sequence we have the following physical considerations: (a) the 
initial fast selective chirp has to be optimally matched to the level structure so that the 
population transfer and selectivity would be extremely high (b) it is necessary to chirp 
up quickly before decay processes become effective and result in false negative counts 
(t c K « 0.16) (c) for t > t c it is necessary to drastically reduce drive strength, since it 
reaches drive strengths which would break the photon blockade through multiphoton 
processes if it persisted. The piecewise linear chirp sequence is fed into a full quantum 
simulation that includes decay, and the 13 parameters of the system and drive are 
optimized with respect to the total readout fidelity. 

The cumulative probability distributions to emit N photons conditioned on starting 
in two initial qubit states are plotted in Fig. |1.8| These distributions were optimized 
for Ti = lfis and then regenerated after varying Ti, in order to depict the effect of 
qubit relaxation on the readout. There are two figures of merit from figure [T~8| one is 
that there exist very high fidelities T — 1 — P(t | t) — P(i | t) exceeding 98% for a 
low threshold around Nth — 20 even for relatively short lived qubits (7\ w 500 ns) . In 
order to take advantage of these fidelities a very low noise amplifier would be needed. 
In addition, we find high contrast and high fidelities (> 90%) for long lived qubits 
Ti > 1.5 /is with thresholds around N t h = 150, which should be accessible with state- 
of-the-art HEMT amplifiers. The limit of obtainable fidelities with this control scheme 
is not due to finite qubit lifetime, as we see from the curve that was simulated for 
Ti = 15 /is. The reason is that after the QCS is generated at t = t c the qubit is 
brought into resonance to form the blockade, and that enhances the anharmonicity 
of the system. Since the drive is not in resonance with all the transition frequencies 
relevant to the wave packet, this leads to the few percent of decay events which are 
unrecovered by the drive. Note also that a useful feature of these distributions is the 
very low level of false positives (red curve for qubit state |t)), for a wide range of 
thresholds, originating from the effectiveness of the photon blockade. 



1.3.3 Experimental considerations 



For experimental applications it is useful to know how robust the fidelity is against de- 
viations of the control pulse parameters from their optimal values. We therefore varied 
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Fig. 1.8 Example photon emission counts histogram generated by quantum trajectories 
simulations for the readout protocol of Fig. |1.7| with optimized parameters. Inset: For a longer 
holding tone, cumulative probability distributions for the number of photons (N) emitted from 
the cavity during the driving time £/, for different qubit decay times (Ti), including a very 
long Ti = lbfis indicating that T\ is not limiting the readout fidelity. For low detection 
thresholds (N th « 20) for distinguishing ft) (N > N th ) from ||) (N < N th ) the fidelity can be 
very high (> 98%) for realistic values of qubit decay in cQED (a few microseconds), and of 
high contrast for more moderate fidelities (> 90%). The distributions also show a almost no 
false positives for higher thresholds Nth > 20 (here tu — 194ns and other parameters as for 
Fig. [T7|). 

the parameters of the initial chirp pulse Sd, c = — 5d, c , $d,q = — 5d, q , £, t c , 
and k independently around their optimal values. Table [Ll] shows for each parameter 
the range of variation for which the fidelity is above 98% (the range cited is the smaller 
of the two ranges above and below the optimal value). The fidelity is most sensitive to 
variations of the duration of the chirp pulse t c which yields a tolerance of ±10% and 
higher ranges (> ±20%) for the rest of the parameters. For achieving slightly less high 
fidelities (> 97%) the bounds for t c increase significantly to ±20%, which is important 
for a realistic experimental setup, since the quench of the Hamiltonian parameters at 
t = t c will take a few nanoseconds with the current microwave technologies. 

The effect of amplifier noise should also be considered. A cryogenic HEMT am- 
plifier with noise temperature of T/v ~ 5if adds noise to the amplitude quadratures 
b x (z,t) = -j=(tf(z,t) ± b(z,i)), b y (z,t) = 4=(tf(z,i) - b(z,t)) of the input signal, 



Parameter 


Optimal value 


Range for T > 98% 


5 dtC (t = 0)/2tt 


-56.0 MHz 


±40% 


5d,c/2K 


21.9 MHz/ns 


±20% 


5 d|g (t = 0)/27r 


-226.2 MHz 


±30% 




-8.2 MHz/ns 


±100% 


£/2tT 


29.4 MHz/ns 


±60% 


to 


10.1 ns 


±10% 


k/2tt 


2.5 MHz 


±60% 



Table 1.1 Values of optimal chirp parameters for achieving maximal fidelity and their rela- 
tive tolerances (given in percents of the optimal value). 



where b(z,t) denotes the annihilation of a photon in the transmission line at posi- 
tion z and time t ( | Clerk et a/., 2 010). The annihilation and creation operators are 
defined such that <fi(zo,t) = (tf (zo,t)b(zo,i)) is the photon flux at the point of en- 
try zq to the amplifier. The dimensionless spectral density of the noise is given by 
S = ksTN/huj « 20 where uo is the frequency of the probe signal. To obtain an opti- 
mal signal to noise ratio (SNR) the signal is measured by time integration, which 
introduces a bandwidth of 1/tf around the carrier frequency, where tf is defined 
in Fig. 1.7 This bandwidth is optimal (Gambetta et a/., 2007) in the sense that it 
lets in all the signal but keeps the noise level minimal (disregarding the short initial 
transient of the chirp). In this model the quadrature noise is assumed to be nor- 
mally distributed with bi y7l0 i se (z,t) ~ A/*(/i = 0, a = \fS~ftf) where i = x,y. It is 
therefore possible to estimate the effect of the amplifier noise by analyzing the dis- 
tribution of the total number of emitted photons, including the photons of noise, 

r t i 2 

i, noise 

quadrature independently. For the cases depicted in Fig. 1.8 where 
SNR « 17 and approximately 3% of errors were added by the noise, which is therefore 
not a significant limit to fidelities in the 90 — 95% range. For T/v > 5 If the effect of 
the noise can be mitigated by increasing tf (up to time 75). 



N = 



2t f l^i 



~-x,y 



where the noise is added to each 
350 we have 



1.3.4 Comparison with other schemes 



For cQED all the necessary components for the above scheme have been experimentally 
demonstrated. Strong qubit-cavity coupling has been demonstrated in many experi- 
ments (iThompson et a/., 19921 iRaimond et a/., 20011 iWallraff et a/., 20041). Strong 



driving of a cavity-qubit system has been shown in (Baur et a/., 2009), with the sys- 



tem behaving in a predictable way, as well as photon blockade (Bishop et a/., 2009) 
and fast dynamical control of the qubit frequency via flux bias lines ( |DiCarlo et ~aL, 
2009). In addition there is evidence both theoretically and experimentally for the in- 
creasing role that quantum coherent control plays in the optimization of these systems 
for tasks of quantum information processing. As examples we can mention improv- 
ing single qubit gates (Motzoi et a/., 2009|, two-qubit gates (Fisher et a/., 2010), and 



population transfer for phase qubits (Jirari et a/., 2009). We therefore believe that the 



readout scheme would be applicable for the transmon, although the control parameters 
would have to be re-optimized due to the effect of additional levels. 



The suggested readout scheme is different from other existing schemes in several 
aspects. Compared to dispersive readout (Blais et a/., 2004) it involves very nonlin- 
ear dynamics and could potentially exhibit much higher fidelity and contrast. Even 
though it relies on a dynamical bistability, it is essentially different from the JBA and 
the scheme from Section 1.3.1 (see also (Bishop et a/., 2010)), since it explicitly oper- 
ates using the quantum photon blockade. Our scheme is also essentially different from 
a recently suggested adaptation of electron-shelving readout to circuit QED ( |Englert| 



et a/., 2010). The latter makes use of a third level in addition to the two levels which 



define the qubit, requires a direct coupling to the qubit with negligible direct driving 
of the resonator, and strong driving of the qubit in the regime where the rotating wave 
approximation breaks down. It is important to stress that the optimization of the con- 
trol parameters in our scheme is only partial, since we have limited ourselves to simple 
linear chirps in this work. Our optimization scheme involves a Simplex algorithm that 
searches for a local maximum for the readout fidelity. Each step the fidelity is calcu- 
lated by solving the stochastic Schrodinger equation for the dynamical evolution for 
different realizations of the measurement signal. More complex modulations (e.g. us- 
ing GRAPE type algorithms) are certainly possible although the standard methods for 
optimal control ( |Khaneja et a/., 20 05) may be difficult to implement here due to the 
large Hilbert space. Therefore we believe that an experimentally based optimization 
using adaptive feedback control (Judson and Rabitz, 1992) might be best option, and 
has the potential to yield superior readout fidelities for higher detection thresholds. 



1.4 Conclusion and future prospects 

In the solid state realization of superconducting qubits, nonlinear oscillator dynamics 
arise naturally from the quantum circuit Hamiltonians, operating in the semiclassical 
or fully quantum regime. Some of the physical phenomena are associated with the 
effect of quantum noise on the oscillator state whereas some are associated with the 
quantized anharmonic ladder of states. In this review we tried to demonstrate the 
theoretical challenges that arise in the high excitation regime of these models and its 
relevance to state-of-the-art experiments. 

Since the experimental methods are evolving rapidly, it becomes intriguing to de- 
velop theoretical schemes of control which utilize larger parts of the accessible Hilbert 
space. This involves extending our understanding of nonlinear response of multiple 
strongly coupled transmon-cavity systems to the high drive power regime. We use an- 
alytical tools and exact simulations to access these regimes and test new protocols. 
It would also be increasingly important to adapt tried and tested optimal control 
techniques to these systems as a support of existing and future experiments. 

This work was supported by the NSF under Grants Nos. DMR-1004406 and DMR- 
0653377, LPS/NSA under ARO Contract No. W911NF-09- 1-0514, and in part by the 
facilities and staff of the Yale University Faculty of Arts and Sciences High Performance 
Computing Center. 
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